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ABSTRACT 

The connection between Gamma-Ray Bursts (GRBs) and their afterglows is currently not well un- 
derstood. Afterglow models of synchrotron emission generated by external shocks in the GRB fireball 
model predict emission detectable in the gamma-ray regime 25 keV). In this paper, we present a 
temporal and spectral analysis of a subset of BATSE GRBs with smooth extended emission tails to 
search for signatures of the "early high-energy afterglow", i.e., afterglow emission that initially begins in 
the gamma-ray phase and subsequently evolves into X-Ray, uv, optical, and radio emission as the blast 
wave is decelerated by the ambient medium. From a sample of 40 GRBs we find that the temporal decays 
are best described with a power-law ~ t^, rather than an exponential, with a mean index (/3) ^ —2. 
Spectral analysis shows that ~ 20% of these events are consistent with fast-cooling synchrotron emission 
for an adiabatic blast wave; three of which are consistent with the blast wave evolution of a jet, with 
Fi, ~ t^'P. This behavior suggests that, in some cases, the emission may originate from a narrow jet, 
possibly consisting of "nuggets" whose angular size are less than l/F, where F is the bulk Lorentz factor. 

Subject headings: gamma rays: bursts 



1. INTRODUCTION 



Afterglow emissions from Gamma-Ray Bursts (GRBs) 
in the X-ray, optical, and radio wavebands are in good 
agreement with afterglow models of relativistic fireballs 



( Wijers, Rees, and Meszaros 1997 . Galama et al. 1998 



Waxman 1997, Vietri 19971). The observed afterglow spec- 



trum is well-described as synchrotron emission that arises 
from the interaction of the relativistic blast wave with bulk 
Lorentz factor Fn ~ 10^ — 10'^ with the ambi ent medium 
( Meszaros and Rees 1997 , Galama et al. 199^ ). The often 
highly variable gamma-ray phase of the burst may reflect 
the physical behavior of the fireball progenitor thr ough col- 
lisions internal to the flow, i.e., int ernal shocks ( Sari an^ 



from internal and external shocks may overlap in time. If 
the internal shocks reflect the activity of the progenitor, 
then the onset of the afterglow may be separated from the 
prompt gamma-ray emission. Since the nature of the pro- 
genitor is not known, the effect of the ambient medium 
on the emission from the progenitor is highly problematic. 
The model therefore does not prohibit internal and ex- 
ternal shock emissions from overlap, while in other cases 
the afterglo w emission may be delaved with respect to the 
GRB (e.g., feari and Piran 1999l|Meszaros and Rees 1999 
Vietri 2000|) 



Piran 1997 



Kobayashi ct al. 1997 ). On the other hand 
Dermer and Mitman (1999) have suggested a blast wave 
with an inhomogenous external medium. Heinz and Begel- 
man (1999) have suggested an inhomogcneous bullet-like 
jet outflow that encounters the interstellar medium 



The detection of optical emission si multaneous with the 
gamma-ray emission of GRB990123 ( Akerlof et al. 1999 ) 
provided the first evidence for two distinct emission com- 
ponents in a GRB; here, the prompt optical emission is 
believed to originate from synchrotron emission in the pro- 
duction of the reverse shock gen erated when the electa en- 
counters the ex t ernal medium (Galama ct al. 199E , |Sar: 



m, ■ 1 i- 1 • V i ii' u 1 /^nn J and Piran 199£, Meszaros and Rees 19991). The gamma- 
Iho procioo roiationohip botwoon the oboorvod GRB and i AwuTuuTToo TT ^ ^ T ^^ r 



the afterglow emission is not well understood. GRBs 
recorded by the BeppoSAX satellite suggest that the X- 
ray afterglow emission may be delayed in time from the 
main GRB (e.g., GRB970228, [Costa et al. 2000D or may 
begin during the GRB emission (e.g., GRB980519, |ln^ 



ray spectrum of GRB990123 can not be extrapolated from 
the spectral flux from the simultaneous optical emission, 
indicating that the optical and gamma ,-rav emission orig- 



inate from two sep arate mechanisms ( Briggs ct al. 1999 



Zand ct al. 199E). In the latter case, it is not clear if the 



Galama et al. 1999|) . Evidence for ov erlapping shock cmis - 
sion was also found in GRB980923 ([Ciblin et al. 1999a|), 



X-ray afterglow is a separate underlying emission compo- 
nent or a continuation of the GRB itself. The internal- 
external shock model presents a scenario in which emission 



where a long power-law decay tail (~ t ^■*) was observed 



in soft gamma-rays (25-300 keV). Two separate emission 
components are favored in this burst because the spec- 
tral characteristics of the tail were markedly different from 
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those of the variable main GRB emission. The spectrum 
in the tail is consistent with that of a slow-cooling syn- 
chrot ron spectnim, similar to the behavior of low-energy 



ergy. 



The resulting synchrotron spectrum of the relativistic 



afterg lows (e.g., Bloom et al. 199^ Vrecswijk ct al. 199£ 

The gamma-rays produced by internal shocks and the 
soft gamma-rays of the "afterglow" may therefore overlap, 
the latter having a signature of power-law decay in the syn- 
chrotron afterglow model. If this is the case, at least some 
GRBs in the BATSE database should show signatures of 
the early external shock emission. These events would 
contain a soft gamma-ray (or hard X-ray) tail component 
that decays as a power-law in their time histories, possi- 
bly superposed upon the variable gamma-ray emission. It 
has been shown that the peak frequency of the initial syn- 
chrotron emission, which depends on the parameters of the 
system (see §2), can peak in hard X-rays or gamma-rays 
( Meszaros and Rees 1992| ). Further, it may be possible to 
see a smoothly decaying GRB that is the result of an exter- 
nal sh ock, i.e., the GRB itself is a "high-energy" afterglow. 



For such GRBs, the subsequent afterglow emission in X- 



rays and optical would then simply be the evolution of the 
burst spectrum. A situation like this might arise when the 
progenitor generates only a single energy release (i.e., no 
internal shocks). 

It is well known that the temporal structures of GRBs 
are very diverse and often contain complex, rapid variabil- 
ity. However, some bursts exhibit smooth decay features 
that persist on timescales as long as, or even longer than, 
the variable emission of the burst. Our investigation fo- 
cuses on the combined temporal and spectral behavior of 
a sample of 40 BATSE GRBs that exhibit smooth decays 
during the later phase of their time histories. Many of 
these events fall into a category of bursts traditionally re- 
ferred to as "FREDs" (Fast Rise, Exponential- like Decay), 
bursts with rapid rise tim es and a smooth extended decay 
dKouvehotou et al. 199"^ ). ]n §2 we present temporal and 
spectral properties of the afterglow synchrotron spectrum. 
In §3 we examine the temporal behavior and spectral char- 
acterisitics of the decay emission for the events in our sam- 
ple and compare their spectra with the model synchrotron 
spectrum. A color-color diagram (CCD) technique is also 
applied to systematically explore the spectral evolution of 
each event. In §4 we present a set of high-energy afterglow 
candidates, followed by a discussion of our results in the 
framework of current fireball models. 

2. SYNCHROTRON SPECTRA FROM EXTERNAL SHOCKS 

Internal shocks are capable of liberating some fraction 
of the total fireball energy Eq = TqMqc^, leaving a sig- 
nificant fraction to be i njected into the extern al medium 
via the external shock (Kobayashi ct al. 1997). However, 
recent simulations sugg est that internal sh ock efficiencies 
can approach ~ 100% (Bcloborodov 2000). Nonetheless, 
as the blast wave sweeps up the external medium, it pro- 
duces a relativistic forward shock and a mildly relativistic 
reverse shock in the opposite direction of the initial flow. 
The reverse shock decelerates the ejecta while the forward 
shock continuously accelerates the electrons into a non- 
thermal distribution of energies described by a power law 
due/d'ye oc 7^7^, where je is the electron Lorentz factor. 
The distribution has a low-energy cutoff given by < 7e- 
Behind the shock, the accelerated electrons and magnetic 
field acquire some fraction, ee and es of the internal en- 



electrons consists of four power-law regions (Sari et al 



1998) defined by three critical frequencies i^a, t'c, and z^m, 
where z^a is the self-absorption frequency, i^c = i^(7c) is the 
cooling frequency, and = i^ijm) is the characteristi c 
synchrotron frequency (see Figure 1 in Sari et al. 1998 ). 
Here, we are only concerned with the high-energy spec- 
trum, therefore we do not consider self-absorption. Elec- 
trons with 7e > 7c cool down to 7c, the Lorentz factor of 
an electron that cools on the hydrodynamic timescale of 
the shock (Piran 1999). The electrons cool rapidly when 
7m > 7c, known as fast-cooling (i.e., > i/^), and cool 
more slowly when 7m < 7c, known as slow- cooling. In the 
fast-cooling regime, the evolution of the shock may range 
from fully radiative (eg ~ 1) to fully adiabatic (ce ^ 1). In 
the slow-cooling mode, the evolution can only be adiabatic, 
since jn-i < 7c- The characteristic synchrotron frequ ency 
of an electron w ith minimum Lorentz factor 7m is ( ^ari| 
and Piran 1999| ) 



= 1.0 X 10 



19 



1/2 



corresponding to a break in the observed spectrum with 
energy 

where F is the bulk Lorentz factor and ni is the constant 
density of the ambient medium. Although the frequency 
in equation 1 depends strongly on the parameters of the 
system, the forward shock may very well peak initially in 
hard X-rays or in gamma rays (Sari and Piran 1999). 

The synchrotron spectrum evolves with time according 
to the hydrodynamic evolution of the shock and the geom- 
etry of the fireball (e.g., spherical or collimated). Specifi- 
cally, the time dependence of Vc and v^i will strongly de- 
pend on the time evolution of the Lorentz factor 7(t). As- 
suming a spherical blast wave and a homogeneous medium, 
for radiative fast-cooling, oc t^^^/^ and i^c oc t^^/^, 
while for adiabatic evolution (fast or slow-cooling) oc 
t~^/^ and Vc The shape of the synchrotron spec- 

trum remains constant with time as Vc and Vrn evolve to 
lower values. In the fast-cooling mode, decays faster 
than i^c, causing a transition in the spectrum from fast to 
slow-cooling. 

Since the break frequencies scale with time as a power- 
law, the spectral energy flux of the synchrotron spectrum, 
Fi, (erg cm~^ keV~^), will also scale as a power-law in 
time so that F^{i', t) oc i'°'t^ , where the spectral and tem- 
poral power-law indices, a and /3, depend on the temporal 
ordering of Vc relative to fm, i-e., fast or slow-cooling. For 
radiative fast-cooling. 



^1/3^-1/3 
^-1/2^-4/7 



-6p)/7 

and for adiabatic fast-cooling, 

r (.1/3^1/6 

F, oc <^ ,.-1/2^-1/4 

I i/-P/2i(2-3p)/4 



(3) 



(4) 
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( Sari et al. 1998 ). For slow-cooling the spectral energy flux 
is 

OC <^ j,-(p-l)/2t-3(p-l)/4 < J. < ly^^ (5) 



(Sari et al. 1998). Note that a simple relation exists be- 
tween the temporal and spectral indices through the value 
of the electron index p for the high-energy spectral slopes 
(v > i>c) and the spectral slope below in the slow-cooling 
regime. Defining the low-energy spectral slope as a and 
the high-energy spectral slope as a' , the following relations 
between the temporal and spectral indices for a spherical 
blast wave are established ( ^ari et al. 1998 ): 



/3 = 



2{6a' - 
3a' /2 
3a/2 
3a72 



l)/7 
-1/2 

-1/2 



(fast-cooling, radiative), 
(fast-cooling, adiabatic), 
(slow-cooling, < ly < Vc), 
(slow-cooling, i^c < ly)- 



(6) 



The numerical value of p is readily determined from the 
measured high-energy spectral slope, p = —2a'. Long 
wavelength afterglow measurements give typical electron 
indices in the range 2.0 < p < 2.5. Although the nature 
of the emission in this model is always synchrotron radi- 
ation with its characteristic slopes and breaks, the time 
dependence of the breaks are affected by the details of the 
geometry and dynamics. 

The relations in equation 6 are only valid in the case 
of a spherical blast wave encountering a constant density 
medium. Rhoads (1999) considered the adiabatic evolu- 
tion of a coUimated or jet-like outflow in which the ejecta 
are confined to a conical volume with a half opening an- 
gle 9c- As the outfiow encounters the external medium, 
the bulk Lorentz factor of the fiow, F, d ecreases with 
radius and time as a power-law (e.g., see Huang et al, 
19990 However, the hydrodynamical evolution of the 
shock changes from a power-law to an exponential r egime 
when 9i, = r-^ ^9c ( |Rhoads 199"9| , ^ari et al. 1999| ). The 
observer is able to discern that the flow is confined to an 
expanding cone rather than a sphere because less radiation 
is observed. In consequence, a break in the light curve to 
a ~ t~P behavior is observed as the ejecta sweep up a 
larger amount of mass. For adiabatic evolution of a jet, 
t'm OC Vc OC t*^ — const , and the peak fl ux scales as 



-Fi^.max OC t-'^ ( [Rhoads 1999| , |Sari et al. 1999| ). Thus the 



spectral flux for an adiabatic jet in the fast-cooling regime 
is given by 



OC 



and for slow-cooling. 



r i.V3i-i/3 

OC <^ i.-(p-i)/2r 

[ u~P/^t~P 



V > Vr . 



(7) 



(8) 



The jet geometry can therefore be tested by the simple re- 
lation (3 = 2a' = —p, irrespective of whether the spectrum 
is fast or slow-cooling, provided that v > max(i'c, i^m)- 

3. ANALYSIS 



We examine the properties of extended decay emis- 
sion in GRBs in the energy range ^ 25-2000 keV us- 
ing data from BATSE, a multi-detector all-sky monitor 
instrument onboard the Compton Gamma-Ray Observa- 
tory (CGRO). BATSE consisted of eight identical detector 
modules placed at the corners of the CGRO in the form of 
an octahedron ( Fishman et al. 198£ ). Each module con- 
tains a Large Area Detector (LAD) composed of a sodium 
iodide crystal scintillator that continuously recorded count 
rates in 1.024 and 2.048 second time intervals with four and 
sixteen energy channels, respectively (known as the DIS- 
CLA and CONT data types). Nominally, a burst trigger 
is declared when the count rates in two or more LADs ex- 
ceed the background count rate by at least 5.5(7. Various 
burst data types are then accumulated, including the four 
channel high time resolution (64 ms) discriminator science 
data (DISCSC). The DISCSC and DISCLA rates cover 
four broad energy channels in the 25-2000 keV range (25- 
50, 50-100, 100-300, > 300 keV). The CONT data span 
roughly the same energy range, but with sixteen energy 
channels and 2.048 s time resolution. 

3.1. Dataset and Background Modeling 

Our dataset was collected by visually selecting events 
from the current BATSE catalog with extended decay fea- 
tures, using DISCSC time histories in the 25-2000 keV 
range. Time histories used in this search had a time res- 
olution of 64 ms or longer, therefore our scan was not 
sensitive to the selection of events from the short class of 



bursts in th e bimodal duration distribution (Kouveliotou 
et al. 1993a). A study of decay emission in short GRBs will 



not be included in this analysis but will be the subject of 
future work. Our search resulted in a sample of 40 bursts, 
17 with a FRED-like profile and 23 that exhibit a period 
of variability followed by a smooth decaying emission tail. 

We grouped events into three categories based on the 
characteristic time history of the bursts: (1) pure FREDs 
(PF) , (2) FREDs with initial variability mainly during the 
peak (FV), and (3) bursts with a period of variability fol- 
lowed by an emission tail (V+T). Note that this catego- 
rization only serves as a descriptive guideline for this anal- 
ysis and does not imply a robust temporal classification 
scheme. Our analysis uses discriminator (DISCSC and 
DISCLA) and continuous (CONT) data from the BATSE 
LADs. 

The source count rates in the ith time bin and the jth 
energy channel, Sij^ were obtained by subtracting the 
background model rates, Bij, from the burst time his- 
tory. The background model rates in the jth energy chan- 
nel were generated by modeling pre and post-burst back- 
ground intervals appropriate for each burst with a poly- 
nomial of order n, where 1 < n < 4. Post-burst intervals 
were chosen at sufficiently late times beyond the tail of the 
burst, since the time when the tail emission disappears into 
the background is somewhat uncertain. This method was 
adequate for bursts with durations less than ~ 200 sec- 
onds. For longer bursts however, the long term variations 
in the background can inhibit knowledge of when the tail 
emission drops below the background level. For this rea- 
son, we applied an orbital background subtraction method 
to events with durations that exceed 200 seconds. This 
technique uses as background the average of the CONT 
data count rates registered when CGRO's orbital position 



4 



is at the point closest in geomagnetic latitude to that at 
the time of the burst on days before and after the burst 
trigger. A complete description of the technique is given 
in Connaughton (2000). 



Contours 



3.2. Temporal Modeling 

In the context of afterglow models, the decay emission 
is usually fit with a power-law. We fit the smooth decay of 
the background subtracted source count rates, Si, in each 
burst with a power-law function of the form 



GRB 970925 
25-300 keV 




RiU) = RoiU - toY 



(9) 



where R{ti) is the model count rate of the ith time bin. 
The free parameters of the model are the amplitude, Rq, in 
counts s~^, the power-law index, /3, and the fiducial point 
of divergence, to, given in seconds. As a general guideline, 
the fit intervals [ti , T2] were selected in a systematic man- 
ner. For the PF bursts, the start time of each fit interval, 
Ti , was taken as the time of half- width at half-maximum 
intensity (HWHM) of the burst. This approach obviously 
does not apply to the V-l-T group of bursts. For these 
events, ri was defined as the bin following the apparent 
end time of the variable emission. The fit interval end 
time, T2, was defined as the time when the amplitude of 
the tail count rates first falls within Icrh^i of the background 
model, where ab,i is the Poisson count rate uncertainty of 
the jth bin of the background model. To avoid obtaining 
a premature value of T2 caused by statistical fluctuations 
in the tail, the amplitude of the count rates in the tail 
was calculated using a moving average of 16 time bins. 
The value of T2 was not particularly sensitive to the width 
of the moving average. We further found that the fitted 
model parameter values were generally insensitive to arbi- 
trarily larger values of T2. 

The model was fit to the data using a Levenberg- 
Marquardt nonlinear least squares minimization algo- 
rithm. The algorithm was modified to incorporate model 
variances rather than data variances in the computation 
of the statistic to avoid overweightin g data points wit h 
strong downward Poisson fluctuations ( Ford et al. 1995| ). 
We performed a set of Monte Carlo simulations to test the 
accuracy of our fitting method. We found a bias in the 
distribution of fitted slopes that was hinged on the cor- 
relation of the {(3, to) model parameters. The bias results 
when the value of ti is too far out in the tail of the power- 
law. In this situation, the curvature of the power-law de- 
cay is undersampled and results in a broad minimum. 
The broad minimum is most easily illustrated by plot- 
ting the joint confidence intervals between P and to. For 
example, Figure 1 shows the Ax^ contours for the fit to 
GRB970925 with Ax^ = 2.3, 6.2, and 11.8, corresponding 
to the 68% (Icr), 95% (2cr), and 99% (3cr) confide nce levels 
for two parameters of interest ( Press et al. 1992] ). 



-10 -5 5 10 

I (seconds) 

Fig. 1. — Ax^ contour plot for GRB970925 that shows the cor- 
relation between /3 and to. Elliptical contours are the Icr (68%), 2a 
(95%), and 3cr (99%) joint confidence intervals. Values of /3 and to 
corrresponding to the fitted minimum are indicated by the filled 
circle. 



Examples of three burst decays from our sample are dis- 
played in Figure 2. The dashed lines indicate the best-fit 
power-law model for each event as listed in Table 1. The 
temporal fit parameters for all events in our sample are 
given in Table 1. The uncertainties in the P and to pa- 
rameters quoted in Table 1 reflect the projection of the 
68% confidence contours onto the axis of the parameter 
of interest, and, in nearly all cases are larger than the 
uncertainties obtained from the covariance matrix in the 
Levenberg-Marquardt algorithm. It is important to point 
out that modeling of the temporal decay of afterglow mea- 
surements at very late times after the burst (e.g., days, 
weeks, and months) in, for example, the optical band does 
not suffer from this bias because the value of to is typically 
set to the trigger time of the burst, i.e., a very good ap- 
proximation to the t rue value of to relativ e to the time of 
the fit interval (e.g., Fruchter et al. 1999). In the case of 
the early decays in GRBs however, the fit is very sensitive 
since we are fitting so close in time to the burst trigger. 

For completeness we also modeled the decay interval 
in each event with an exponential function of the form 
Re{ti) ~ exp[— (ti — Ti)/re] with the amplitude and expo- 
nential decay constant Tg as free parameters. We find that 
only 12 of the 41 fits resulted in a lower reduced x^ value, 
Xr, than that of the power- law model. The largest value 
of Axr of these events was only 1.15 while most other 
events had Axr ~ 0.3 or less, indicating that the power- 
law is nearly as good a fit as the exponential. For events in 
which the exponential model was a poor fit, the power-law 
fits were strongly favored with Axr values as high as 7.5. 
Our results are consistent with that of a similar study for 
a small number of GRBs perform ed early in the BATSE 
mission ( ^chaefer and Dyson 1995 ). 
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Table 1 

Summary of Temporal Fits (25-300 keV) 



GRB 


Trigger 


Time Profile^i) 


[ti,T2] (s) 


/3 


to 


d.o.f. 


X^/d.o.f. 


910602 


257 


FV 


[19,361] 


— 1 74 


— 13 69 •^•^■00 


166 


1.41 


910814c 


676 


V+T 


[67,109] 


„ 0.65 
+[).4() 


58.05^2 ()5 


39 


0.98 


910927 


829 


FV 


[15,52] 




7 69^2-" 

' -^-^ + 2 05 


34 


1.20 


911016 


907 


FV 


[140,292] 


< -2.06 


< 113.73 


147 


4.03 


920218 


1419 


FV 


[130,204] 


-2 1^-0-54 
^•^^+0.33 


120.90;i;02 


71 


1.32 


920502 


1578 


FV 


[14,99] 


—9 ^1 —3-27 




81 


1.04 


920622 


1663 


V+T 


[21,80] 


< -2.98 


< 4.77 


56 


1.75 


920801 


1733 


PF 


[9,123] 


— 1 91 


9 AO 

—9 S4 

^•°^+3.06 


110 


2.15 


920813 


1807 


FV 


[39,150] 


—2 in-6.05 


18.07;Ji|7i 


106 


1.57 


920901 


1885 


PF 


[19,200] 


< -2.02 


< -12.81 


175 


1.15 


yzizu / 




r V 


[iZ,ODj 


o as-0.03 
-2.»8^0.02 


74-0.03 
' +0.02 




1 -3 


yoUlUo 




r V 




1 ^^"O-SS 
-l-'O+0.12 


1 OQ — 9.23 

-12.28_^g Og 


01 Q 


n QQ 
U.oo 






V + i 


[3,96] 


n 71 -0.31 


9 .1 -4.76 
^•^^+2.12 


89 


0.91 


930612 


2387 


PF 


[17,264] 


-2.09+»:f; 


T _Q 46 
-^-^-^ + 1.48 


190 


1.24 


931223 


2706 


PF 


[10,72] 


_o 1:7-6.71 
^•'-''+0.10 


-4 7=5-60.4 
^-'^+1.03 


59 


1.26 


940218 


2833 


FV 


[7,50] 


n ns 
fjQ — u.uo 

■'•^°+a.04 


^-"■'+0.03 


40 


1.60 


940419b 


2939 


PF 


[17,264] 


< -1.75 


< —17.48 


239 


1.51 




3257 


OTP 


[17,131] 


< —1.95 


< —8.89 


109 


1.67 


OK1 1 OA 


ooyo 


V + 1 


[no QQl 


1 07-0.22 
1+0.02 


14 40-0.05 
^^•^^+1.03 


1 


l.Uo 


960530(1) 


5478 


PF 


[9,121] 


J--^y-|-0.04 


2.36_|_Q;2g 


108 


0.88 




5478 


PF 


[273,515] 


< —2.13 


< 251.46 


234 


1.33 


970302 


6111 


PF 


[9,133] 


-i-49+o:i'-i 


„ .... 1 n 4Q 

-o-8i+;:'if 


119 


1.85 


970411 


6168 


FV 


[44,398] 




16-01+423 


171 


0.95 


970925 


6397 


PF 


[12,85] 




-i-65;^:25 


69 


1.71 


971127 


6504 


PF 


[11,198] 


-1 96""*'^ 


-4 26"^-*'^ 
^-^"+1.57 


180 


1.44 


971208 


6526 


Fb 


[361,2995] 


1 Q^-0.11 
-1-34+0.01 


-26.34 55 

27.65;^:?; 


1047 


2.32 


980301 


6621 


PF 


[35,87] 


^•^"+0.05 


48 


1.37 


980306 


6629 


FV 


[242,402] 


-1 71-0-71 

^-'^+0.50 


219.03;5-« 


154 


1.09 


980325 


6657 


PF 


[20,152] 
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^-•^^+0.04 
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'-^^+1.58 
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2.36 


980329 


6665 


V+T 


[16,59] 


< -2.39 


< 6.64 


40 


2.48 


981203 
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FV 


[60,851] 


-1 61-0-"is 
-^-"-^+0.002 


_17 1-6.15 
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i.oa^]^ 20 
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-2 14^1-^''' 


7-15.71 

' +2.0 
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1.14 


990220 
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PF 


[17,196] 


_1 07-1.04 

-^-=^'+0.30 


n '70 — 16.86 
^y-'^'+fi.lO 


85 


0.83 


990316 
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PF 


[25,158] 


< -2.20 


< 11.06 


128 


1.54 


990322 


7488 


PF 


[4,110] 


-U.S 1^+0.01 


-, 97-0.05 

l-^'+0.13 


102 


1.17 


990415 


7520 


V+T 


[44,122] 


< -2.23 


< 29.44 


75 


1.54 


990518 
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V+T 
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1.61 
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_9 qf; — 1.15 
^•■-''^+0.03 
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(1) Abbreviations for time profile descriptions: PF = pure FRED (i.e., single smooth pulse), FV = FRED 
-with multiple pulses near the peak, V+T = variability + tail. 

(2) First emission episode of trigger 5478. 
Second emission episode of trigger 5478. 

(4) The longest duration GRB observed by BATSE (~ 3000 seconds). 
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Fig. 2. — A logarithmic plot of the time histories of three events 
in the 25-300 keV range. The dashed line is the best-fit power-law 
model for each burst. The time intervals for each fit are listed in 
the fourth column in Table 1. 



3.3. Spectral Modeling 

Nearly all GRB spectra are adequately modeled with a 
low and high-energy power-law function smoothly joined 
over s ome energy range within the BATSE energy band- 



pass (Band et al. 1993, Prcccc ct al. 200C). Curvature 



in the spectrum is almost always observed, although on 
rare occasions a broken po wer-law (BPL) mo del is a bet- 
ter represention of the data ( |Preece et al. 199^ ). Spectra of 
X-ray afterglows in the 2-10 keV range observed with Bep- 
poSAX are best fit with a single powe r-law, with spectra l 
indices that range from -1.5 to -2.3 ( [Costa et al. 20001 ). 
Recently it has been noted that the breaks in the syn- 
chrotron spectrum may not be sharp, but rather smooth 
(Granot and Sari 2001). Therefore we chose two spectral 
forms to model the spectra of the gamma-ray tails: a single 
power-law as a baseline function, and a smoothly broken 
power-law (SBPL). The SBPL was chosen to enable di- 
rect comparison with the spectral form of the synchrotron 
shock model. 

Because we are interested in the spectral behavior dur- 
ing late times of the burst, the CONT data-type from in- 
dividual detectors is the optimum choice of available data- 
types from the LADs. CONT affords the best compromise 
between temporal and energy coverage, with 16 energy 
channels and 2.048 s time resolution. Coarse temporal 
and energy bins are required as we are dealing with a sig- 
nal that continuously decays with time. 

We model the photon spectrum (photons s^^ cm^ 



keV~^) using the standard deconvolution and the 
Levenberg-Marquardt nonlinear least-squares fitting algo- 
rithm that incorporates model variances. The spectra were 
modeled using CONT channels 2-14, which covered the en- 
ergy range of ~ 30-1800 keV. Count spectra from the two 
brightest detectors (i.e., the two detectors with the small- 
est source angles to the LAD normal vector) were generally 
used to make the fits. In some cases, the source angles dif- 
fered substantially (> 20°) resulting in a normalization 
offset in the count spectra between the two detectors. The 
data from the two detectors were fit jointly with a mul- 
tiplicative effective area correction term in the spectral 
model. For bursts in which the effective area correction 
was small (^ 5%), we summed the CONT count rates from 
the individual detectors to maximize the count statistics. 
These were cases in which the source angles of the two 
detectors differed by only a few degrees. Detectors with 
angles to the source exceeding 60° or with strong signal 
from sources such as Vela X-1 or Cyg X-1 were excluded 
from the fit. 

The free parameters of the power-law spectral model 
are the amplitude and the power-law index ap. The free 
parameters of the BPL and SBPL are the amplitude, low- 
energy index aiow, high-energy index ahigh, and the break 
energy Eh. The slopes of the spectral energy flux, are 
readily obtained from the simple relations: a — aiow + 1 
and a' — ahigh + 1- 

For the decay emission of each burst, we modeled the 
time-integrated spectrum defined over a time interval that 
was either the same as or shorter in length than the time 
interval used in making the temporal fits. In all cases, 
the time interval was restricted to the region of the burst 
during the power-law decay. Shorter time intervals were 
used for events with weaker signal-to-noise. For the PF 
class bursts, however, we selected the entire burst emis- 
sion (spectra of the burst intervals starting at the peak 
gave nearly identical parameter values). For the majority 
of events, the SBPL model was preferred over the single 
power-law model. Summarized in Table 2 are the spectral 
fit parameters for events where the SBPL was the better 
choice of model based on the statistic. In 6 bursts, 
however, the fitted value of the high-energy slope was un- 
usually steep (^ —4.0). Preece et al. (1998) pointed out 
that spectral models with curvature can sometimes over- 
estimate the steepness of the high-energy slope, depend- 
in g on how well the data tolerate curvature (see figure 1 



Preece et al. 1998|). The broken power-law model was 



therefore used for these 6 events, resulting in slightly bet- 
ter reduced values and better constrained values of the 
fitted high-energy slope. This resulted in spectral param- 
eters for a total of 20 bursts. The reduced x^ values are 
reasonable, although a few bursts for which joint fits were 
made tended to give slightly larger values {x^/d.o.f. ^ 2). 
Given in the table are the fitted values of the spectral in- 
dices, the break energy, and their uncertainties from the 
covariance matrix. Also given is the difference in spectral 
slope across the break energy, A =\ a' — a \ , and the value 
of p calculated from the high-energy spectral slope. 

Of the remaining bursts, 9 events were best represented 
by the single power-law function. The best fit param- 
eters and the corresponding value of p are presented in 
Table 3. The spectral fits for the remaining 12 events re- 
sulted in poor x^ values and poorly constrained parameter 
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Table 2 

Time-Integrated Spectral Fits: Smoothly Broken Power-Law 



GRB 


a(2) 


a' (2) 


Eb (keV) 


x'^/d.o.f. 


A 


P 


category^*) 


910927 


0.06 ± 0.04 


-3.27 ±0.17 


158 ±7 


1.84 


3.33 ±0.17 


6.54 ±0.17 




920218 


-0.69 ± 0.01 


-1.45 ±0.05 


175 ± 7 


2.37 


0.76 ± 0.15 


2.90 ± 0.05 




920502 


-0.10 ± 0.05 


-1.96 ± 0.20 


183 ± 20 


1.22 


1.86 ± 0.21 


3.92 ± 0.20 




920622 


-0.49 ± 0.04 


-1.48 ± 0.26 


263 ± 68 


1.97 


0.99 ± 0.26 


2.96 ± 0.26 


(i).(iii) 


920801 


-0.23 ± 0.05 


-0.98 ± 0.15 


252 ± 70 


1.13 


0.75 ± 0.16 


1.96 ± 0.15 


(iii) 


930612 


-0.01 ± 0.17 


— 1.46 ± 0.12 


108 ± 19 


0.49 


1.45 ± 0.21 


2.92 ± 0.12 


931223 


-0.24 ± 0.06 


-1.39 ±0.13 


141 ± 12 


1.28 


1.15 ±0.14 


2.79 ±0.13 




940419b 


-0.52 ±0.17 


-1.66 ±0.55 


140 ± 73 


1.15 


1.14 ±0.58 


3.32 ± 0.55 


(i) 


941026 


-0.24 ± 0.03 


-1.55 ±0.08 


150 ±6 


2.32 


1.32 ± 0.09 


3.10 ±0.08 




960530(1) 


-0.10 ±0.18 


-1.48 ±0.22 


128 ± 33 


1.59 


1.47 ±0.28 


2.96 ± 0.22 




960530(2) 


-0.56 ± 0.04 


-1.88 ±0.37 


203 ± 25 


1.54 


1.33 ± 0.37 


3.76 ± 0.37 


(i)(5) 


970411 


-0.35 ±0.04 


-1.10 ±0.08 


228 ± 41 


1.60 


0.75 ± 0.09 


2.20 ±0.08 


(iii) 


970925 


-0.02 ±0.32 


-1.28 ±0.13 


101 ± 32 


1.10 


1.26 ± 0.35 


2.56 ± 0.13 


(iii) 


971208 


-0.55 ±0.02 


-2.03 ±0.04 


179 ±6 


2.02 


1.48 ± 0.04 


4.06 ± 0.04 


980301 


-0.55 ±0.22 


-1.24 ±0.13 


76 ± 34 


0.80 


0.69 ± 0.13 


2.48 ±0.13 


(i),(iii) 


981203 


0.29 ± 0.05 


-0.59 ±0.01 


124 ±8 


2.28 


0.88 ± 0.05 


... (3) 


(iii) 


990102 


0.53 ± 0.16 


-1.82 ±0.15 


121 ± 14 


1.03 


2.35 ± 0.22 


3.64 ±0.15 


990220 


0.64 ±0.24 


-1.36 ±0.08 


94± 13 


5.79 


2.00 ± 0.25 


2.72 ± 0.08 




990316 


-0.58 ±0.04 


-1.52 ±0.09 


145 ± 18 


3.44 


0.94 ±0.10 


3.04 ±0.09 




990518 


-0.52 ± 0.04 


-1.45 ±0.20 


174 ± 22 


2.07 


0.93 ± 0.20 


2.90 ± 0.20 


(i) 



('^^For most events the fluence interval is the same interval used in making the temporal fit. 

(^^Here, a and a' are the indices of the of the spectral energy flux, (erg s~^ cm"^ keV~^), i.e. a = aiow + 1 and a' = ahigh + 1- 

(^^Note that the value of a is consistent with the spectral slope below Uc in the fast-cooling mode and the spectral slope below Um in the 
slow-cooling mode. In the former case, p is undetermined. For slow-cooling, p = 2.18. 

Characterisitc signatures of the synchrotron spectrum as described in §3.3 of the text. 

(5) Within 1.5 sigma 

Note. — Uncertainties in the model parameters axe taken from the covariance matrix. 



Table 3 

Time-Integrated Spectral Fits: Single Power-Law^^^ 



GRB 


Trigger 


ap(2) 


xVd.o.f. 


p(3) 


910602 


257 


-0.52 ± 0.01 


1.05 


1.04 ±0.01 


911016 


907 


-1.77 ±0.22 


1.22 


3.54 ±0.22 


920813 


1807 


-1.01 ± 0.02 


13.2 


2.02 ± 0.02 


930131 


2151 


-0.89 ±0.11 


0.89 


1.78 ±0.11 


951104 


3893 


-1.47 ±0.05 


4.96 


2.94 ± 0.05 


970302 


6111 


-0.73 ±0.10 


1.13 


1.46 ±0.10 


990322 


7488 


-0.42 ± 0.06 


1.51 


0.84 ± 0.06 


990415 


7520 


-1.28 ±0.16 


1.19 


2.56 ±0.16 


000103 


7932 


-1.99 ± 0.43 


1.34 


3.98 ± 0.43 



(1) por most events the flucncc interval is the same interval used in making the temporal fit. 

(2) Powcr-law index is the index of the spectral energy flux, Fi, (erg s~^ cm"'^ keV~^). 

(3) Here, p = -2ap. 



Note. — Uncertainties in the model parameters are taken from the covaxiance matrix. 
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values, regardless of the choice of spectral model. These 
are clearly cases when the counting statistics are too poor 
to constrain the model parameters and therefore were ex- 
cluded. The results in Table 3 should be interpreted with 
some degree of caution. These events may be cases in 
which the flux level was to low, causing the break in the 
spectrum to be washed out in the counting noise. In such 
cases, the single power-law will often be adequate to model 
the spectrum, even though the true burst spectrum may 
contain a break. 

A careful inspection of Table 2 immediately allows 
us to identify high-energy afterglow candidates based on 
three characteristic signatures of the synchrotron spec- 
trum which we categorize as the following: (i) in the fast- 
cooling mode, the spectral slope below the high-energy 
break (vm) is always —1/2 for radiative or adiabatic evo- 
lution, as seen from equations 3 and 4, (ii) in the slow- 
cooling mode, the change in spectral slope across the high- 
energy break (vc) is always 1/2, as seen from equation 5, 
and (iii) the electron energy index, p, calculated from the 
measured spectral slope should have a value in the range 
2.0 < p < 2.5, the typical range derived from afterglows 
observed at X-ray, optical, and radio wavelengths. 

Applying these criteria, we label events with these prop- 
erties in the last column of Table 2. We thus immedi- 
ately identify several fast-cooling candidates: GRB920622, 
GRB940419b, GRB960530(2), GRB980301. Each of 
these events has a value of a within one-sigma of —0.5 
and a value of p similar to those found for afterglows. 
GRB970411, GRB971208, GRB990316, and GRB990518 
are only marginally consistent with fast-cooling, having 
larger p values and reduced x^'s. None of the events in 
Table 2 are consistent (within one-sigma) with A = 0.5, 
suggesting no slow-cooling candidates (however, 5 events 
[GRB920218, GRB920622, GRB920801, GRB940419b, 
and GRB980301] have values within two-sigma). A total 
of 9 bursts in Table 2 are ruled out as high-energy after- 
glow candidates because their spectral parameters bear no 
resemblance to the fast or slow-cooling synchrotron spec- 
trum. One event of notable interest is GRB981203, which 
has a and a' values consistent with a cooling break Vc , as 
opposed to Vni, in the fast-cooling spectrum. This implies 
a i/yn break above ~ 2 MeV, while the value of p remains 
unconstrained by the data. In §4 the early high-energy 
afterglow candidates are discussed in greater detail. 

Obviously, for the single power-law events listed in Table 
3 we have less spectral information. The value oip given in 
Table 3 is derived from ap under the assumption that ap is 
the slope above the break, for fast or slow-cooling. Clearly 
this need not be the case. A case in point is GRB910602, 
which has ap ~ —0.52 ± 0.01, a value consistent with the 
spectral slope of fast-cooling for i/c < < i^m- In this 
interpretation the cooling break, would be below the 
BATSE window and above. The value of p would be 
undetermined from the data. Scanning the values of p 
given in Table 3, we find p = 2.56 ± 0.16 for GRB990415, 
a typical value for afterglows. This suggests that the mea- 



expected value if p = 2.5. 



sured slope ap 



-1.28 ± 0.16 could be the slow or fast- 



cooling high-energy slope. What values of the spectral 
slope do we expect to observe below Vc for slow-cooling? 
If we assume a value of p = 2.5, then the calculated slope 
below I'c is a = —{p — l)/2 = —0.75. We find one burst, 
GRB970302, with a„ = -0.73 ± 0.10, consistent with the 



3=3a72+l/2 
(S = 2{6a'+l)/7 
S = 2a' = -p 
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Fig. 3. — A plot of high-energy spectral index vs. temporal index 
for the twenty events listed in Table 2. Also plotted are the linear 
relationships expected from the evolution of an adiabatic spherical 
blast wave {solid line, dash-dot line) and a jet {dashed line). The 
five bursts labeled on the plot {fill diamonds) are within one-sigma 
of the f3 = 2a' line. 

An additional constraint we can apply to the data is a 
comparison of the measured temporal slopes with their ex- 
pected values derived from the measured spectral indices 
given in the expressions in equation 6. A plot of tem- 
poral vs. spectral index for the data in Table 2 is shown 
in Figure 3. Here, the spectral index is the high-energy 
spectral energy index, a', in the third column of Table 
2. For comparison with the models, we plot the possible 
linear relationships between a' and /3 given in equation 6. 
Note that this plot should be interpreted with a certain 
degree of caution. The expected values of a' are some- 
what restricted by the possible range of p values between 
2.0 and 2.5 predicte d by Fermi accelerati on models (e.g.. 
Gallant et al. 1999], [Gallant et al. 20001). Interestingly, 



however, five events {filled diamonds) are consistent with 
the P = —p line for adiabatic jet evolution. We address 
this implication in detail in §4 and §5. 

A similar plot is shown in Figure 4 for the single power- 
law fits from Table 3. In general, all but one event ap- 
pear consistent with the relations between the temporal 
and spectral indices expected from external shocks. Thus, 
closer inspection of the spectral parameters (e.g., ap and 
p) in Table 3 is required to establish if these are viable 
high-energy afterglow candidates (see §4). 

3.4. Spectral Evolution: Color-Color Diagrams 

The evolution of the synchrotron spectrum is unique 
for a given hydrodynamical evolution of the blast wave. 
This evolution can be traced in a graphical form using 
a Color-Color Diagram (CCD). The CCD method is a 
model-independent technique that characterizes the spec- 
tral evolution of the burst over a specified energy range. 
With this method, a comparison of spectral evolution pat- 
terns among GRBs can be made in addition to a com- 
parison with patterns expected from the evolution of the 
synchrotron spectrum. 



-1 




The crescent patt e rn is typical among GRBs ( Kouve 
liotou et al. 1993b|, piblin et al. 1999b| ), however the 
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Fig. 4. — A plot of spectral index vs. temporal index for the 
9 events listed in Table 3. The asymmetric uncertainties in the 
temporal indices reflect the 68% joint confindence intervals between 
f3 and to- Thus the uncertainties for bursts with highly elongated 
Ax^ contours are indicated by downward arrows. Lines indicate the 
linear relationships expected in the synchrotron afterglow spectrum: 
i^m < u < Uc in slow cooling, (3 = 3ap/2 {dotted line), Vrn < v in adi- 
abatic fast-cooling or Uc < v in slow cooling, f) = 3ap/2-|- 1/2 (solid 
line), and for Um < u in radiative fast cooling, /3 = 2(6«p + l)/7 
(dash-dot). Also plotted is the (3 = 2ap line for adiabatic jet evo- 
lution. Since we do not observe a break in the spectrum for these 
events and thus can not distinguish if the observed power-law index 
is the low or high-energy index, we plot all possible relationships 
between the spectral and temporal index. 



The CCD is a plot of the hard color vs. the soft color, 
where the hard and soft colors are defined as the hard- 
ness ratios (i.e., ratios of the count rates) between (100- 
300 keV/50-100 keV) and (50-100 keV/25-50 keV), respec- 
tively. To construct the CCDs we use the count rates in 
the three lowest (25-300 keV) of the four broad energy 
channels from the DISCSC data. We select a time inter- 
val large enough to cover most of the burst emission until 
the statistical noise begins to dominate. These bins are 
identified by hardness ratios with two-sigma upper limits. 
We also fold the fast and slow-cooling broken power-law 
synchrotron spectra of a spherical blast wave through the 
LAD detector response to obtain the expected count spec- 
trum (shown as the dashed and solid lines, respectively, 
in Figures 5-8). We assume the fast-coohng spectrum is 
radiative, with p = 2.4, aiow = —1.5, and allow E^n to 
evolve from 220 25 keV. For the slow-cooling spectrum, 
we also assume p — 2.4 and the same evolution for E^. 

Giblin et al. (1999a) have shown that the tail emission 
from GRB980923 resembles that of afterglow synchrotron 
emission due to an external shock. To illustrate the useful- 
ness of the CCD technique, we show in Figure 5 the CCD 
for GRB980923. In this representation, the time evolution 
of the burst is preserved by a color sequence of the hard- 
ness ratios, with black/violet/blue signaling the onset of 
the burst and yellow/red signaling the end of the burst. 
The left panel shows the CCD for the time interval that 
brackets the entire burst (variability -I- tail). The variable 
emission of the burst shows a crescent-like pattern decou- 
pled from a cluster of points that represent the tail of the 
burst. 



clustering is less common. The crescent track exhibits a 
sawtoothing of soft-hard-soft evolution, indicative of the 
spectral behavior of the individual pulses that comprise 
the main burst emission. The pattern drastically changes 
when the variability ceases and the tail becomes visible. 
The tail cluster overlaps the region of the two-color plane 
that contains the evolution of the slow and fast-cooling 
synchrotron spectrum. This is best illustrated in the right 
panel of Figure 5, where the CCD is constructed from a 
longer time interval in the tail only. Unfortunately the 
CCD pattern of the tail is not completely resolved due to 
the increasingly large uncertainties in the hardness ratios 
that arise from the decreasing flux level. However, the 
points do lie in the correct region of the diagram. This 
decoupling of the points in the model-independent CCD 
is clear evidence for two distinct spectral components ob- 
served in a GRB. 

Figure 6 shows the CCDs for the four fast-cooling can- 
didates identified based on their spectral parameters in 
Table 2. The pattern for GRB920622 bears a striking re- 
semblance to that of GRB980923 in the left panel of Fig- 
ure 5. Like GRB980923, this burst contains a period of 
variability followed by a very smooth emission tail. The 
crescent pattern that arises from the variable part of the 
burst is clearly visible and spans nearly the same range of 
soft and hard color indices. However, the clear disconti- 
nuity between the burst emission and the tail emission in 
the CCD of GRB980923 is not as pronounced in the CCD 
of GRB920622. Nonetheless, the tail emission (denoted by 
the orange and red points) lies in the same region as those 
of GRB980923 and the synchrotron afterglow spectrum. 
The CCD pattern for GRB949419b (in the PF class) is 
similar, but appears more cluster-like in the region of syn- 
chrotron evolution. 

Another burst of interest in the PF category is 
GRB960530. This event consists of two FREDs separated 
by ~ 200 seconds. The 2nd FRED only has about half 
of the peak intensity as the first. The CCD for the 2nd 
FRED is seen in Figure 6. The episode begins very hard 
on the rise and evolves through the synchrotron spectrum 
during the decay. Interestingly, although the 1st episode of 
GRB960530 is also a FRED, its color-color diagram (Fig- 
ure 8) shows a broad crescent pattern that evolves much 
farther away from the synchrotron pattern. This may be a 
case where the external shock is clearly decoupled in time 
from the GRB. 

The last event in Figure 6, GRB980301, shows an in- 
triguing pattern that closely resembles the evolution pat- 
tern of the synchrotron spectrum. The evolution is mainly 
shaped like a reverse "L" , but marginally offset to higher 
soft color values and lower hard color values. In this case 
it is difficult to argue in favor or against the synchrotron 
model. 

Figure 7 shows the CCDs of the four events from Table 
3 that are only marginally consistent with fast-cooling. 
GRB970411 and GRB990518 show similar patterns that 
resemble those in Figure 6, but the consistency with the 
synchrotron pattern is weak. GRB971208 shows little evo- 
lution and a cluster of points partially overlapping the 
synchrotron region. GRB990316 shows a nearly identical 
pattern to that of GRB980301. 
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Fig. 5. — Left panel: Color-color diagram of GRB980923 variability episode and tail emission (0-70 seconds). The crescent-like pattern is 
the variable emission while the decoupled cluster of points is the tail emission. Right panel: Color-color diagram of the tail emission only, 
covering a much longer time interval (40-100 seconds). Also plotted are the evolution patterns expected from the slow (solid) and fast-cooling 
(dashed) synchrotron spectrum. 
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Fig. 6. — Color-color diagrams for the four fast-cooling candidates from Table 3. Arrows indicate two-sigma upper limits. 
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Fig. 7. — Color-color diagrams for the marginal fast-cooling candidates from Table 3. Arrows indicate two-sigma upper limits. 



A total of 23 events from our sample showed CCDs in- 
consistent with the evolution of the synchrotron spectrum. 
Rather, most of these events showed the crescent pattern 
that are common among GRBs, as depicted in Figure 8. 
These events span a much larger range of hard and soft 
colors than expected from the synchrotron emission alone. 
Others were too weak to distinguish a pattern. 

4. HIGH-ENERGY AFTERGLOW CANDIDATES 

We identify a total of 8 from our sample of 40 events as 
high-energy afterglow candidates based on their observed 
spectral parameters and color-color diagrams. Each burst 
is discussed in detail below. 

4.0.1. GRB 910602 

The observed spectral slope for GRB9 10602 is consis- 
tent with the spectral slope below in the fast-cooling 
spectrum, although the slope below in slow-cooling can 
not be ruled out, implying a value of p = 2. A series of 
time-resolved fits with a uniform time resolution of 4.096 
s revealed no softening of the spectrum, i.e, the slope re- 
mained constant with ap ~ —0.5 throughout the tail. Ap- 
plying the relations in equation 6, for slow-cooling we ex- 
pect (3 = —0.78 ±0.01 and, from equations 4 and 5, for 
fast-cooling we expect (3 — —4/7 (radiative) or /3 ~ —1/4 
(adiabatic). These values do not agree with the measured 
value (3 = — 1.74^q;J^. Although the spectrum appears 
to be consistent with that of the synchrotron spectrum, 
the evolution does not appear to be consistent with the 
evolution of a spherical blast wave. 



4.0.2. GRB 920622 



The time history of this burst bears a striking resem- 
blance to that of GRB980923 reported by Gibhn et al. 
(1999). Initially the burst is highly variable, then at ^ 18 
s after the trigger time the burst enters a phase of smooth 
decay that lasts until ^ 50 s after the trigger. From Ta- 
ble 2, the time-integrated spectral fit suggests fast-cooling, 
with low-energy index a = —0.49 ± 0.04. From the high- 
energy index, a value of p = 2.96 ± 0.26 is inferred. The 
value of A = 0.99 ± 0.26 is marginally consistent (within 
two-sigma) with the expected value of 0.5 for slow-cooling. 
The time-integrated spectrum of the variable emission of 
the burst is in contrast with the fluence spectrum of the 
tail. The variable emission gives = —0.07 ± 0.01, 
a; = -1.50±0.04, and E^,b = 370±12 keV, suggestive of a 
spectral change near ~ 18 seconds. Note that the spectral 
parameters of the variable emission are not consistent with 
the synchrotron spectrum. We binned the tail emission 
into three time bins with S/N > 45 to model the spectral 
evolution, however the parameters were poorly constrained 
due to the steep nature of the flux decay. The spectral evo- 
lution of the variability -f tail emission, however, can be 
seen in Figure 6. The tail of the burst appears consistent 
with the region of the diagram defined by the evolution 
of the synchrotron spectrum. The measured temporal in- 
dex of the tail is (3 < —2.98, clearly inconsistent with the 
expected values for (3 for v < v^. Interestingly, however, 
P is nearly identical to the value of p inferred from the 
high-energy slope, as expected for jet evolution. 
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Fig. 8. — A sample of color-color diagrams of GRBs that are not consistent with the evolution of the synchrotron spectrum. Arrows indicate 
two-sigma upper limits. 



4.0.3. GRB 940419b 

The smooth rise and decay structure of this burst 
place it in the PF category. Like the tail emission of 
GRB920622, this burst also shows a low-energy slope con- 
sistent with the fast-cooling spectral slope below ly^i- The 
measured temporal slope is (3 < —1.75. While this slope 
is not consistent with the temporal index below i^^, it is 
consistent with the expected values of /? = —2.56 ± 0.55 
(radiative) and /3 = — 1.95 ± 0.55 (adiabatic) for > i/m, 
given the large uncertainties. We further binned the data 
in the tail to S/N > 45 and constrained the evolution 
of the break energy by holding the low and high-energy 
spectral indices fixed to their values derived from the time- 
integrated fit. We find Eb - {t-to)-^-^^^°-^^ for to fixed at 
-25.0 s (x^/d.o.f. = 7.52/8). Additional fits with other 
larger values of to gave slightly shallower indices, as ex- 
pected if one compares the behavior to that of f3 and to- 
For an adiabatic fast-cooling spherical blast-wave we ex- 
pect the break energy to decay as —1.5, within two-sigma 
of our measured value. As seen from the CCD in Fig- 
ure 5, the spectral evolution of this burst is very close to 
that of the synchrotron spectrum, although the rise of the 
burst tends to be somewhat harder in the soft color index 
than expected from evolution of the synchrotron spectrum 
alone. 

4.0.4. GRB 960530 

GRB960530 is of particular interest because of its strik- 
ing temporal behavior. The burst has two distinct episodes 
of emission, each having a FRED-like time profile. The 



second episode, much weaker with a peak intensity less 
than half of the first, occurs ~ 200 seconds after the first. 
As seen in Table 2, the low-energy slope of the second 
episode is consistent with the fast-cooling slope below fm, 
however the value of p — 3.76 ± 0.37 is large mainly be- 
cause the high-energy slope is not well-constrained. The 
value of p derived for the first episode is not unreasonable, 
however the low-energy slope is roughly three-sigma away 
from the expected value of —0.5. The decay index for the 
second episode is /3 < —2.13, not inconsistent with the 
expected values /3 = -4.57 ± 1.57 and (3 = -3.75 ± 1.57 
for radiative and adiabatic fast-cooling, respectively. The 
CCD of the second emission episode of this burst (Figure 
6) indicates that during the decay the emission evolves 
into the synchrotron spectrum. 

4.0.5. GRB 970411 

From Table 2, the low-energy index of this burst is 
nearly four-sigma from the value —0.5 expected in the 
fast-cooling regime. However, it does have p = 2.2 ± 0.08, 
consistent with typical afterglow values and particle acceL 



eration models of relativistic shocks (e.g.. Gallant et al 



2000). Additionally, the change in slope, A is less than 



three-sigma from the expected value of —0.5 for the cool- 
ing break in the slow-cooling regime. For slow-cooling, 
we expect the temporal slope to be /3 = —0.53 ± 0.04 for 
V < Vc and fi = —1.15 ± 0.08 for v > v^- For v > in 
fast-cooling we expect (3 = —1.60 ± 0.08 (radiative) and 
P = —1.15 ± 0.08 (adiabatic). Our measured value of the 



decay, /3 



-2.06 



-0.21 
-1-0.20' 



is marginally consistent (within 
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two-sigma) with the radiative fast-cooUng slope. More no- 
tably, it is consistent (within one-sigma) with the value 
of p. A series of spectral fits during the tail of the burst 
holding the low and high-energy spectral indices constant 
show that Eh decays with time described by a power- 
law of the form Eb ^ {t - to)-0-96±o.26 ^j. to ^ 16 s 
(x^ /d.o.f. = 4.02/3), marginally consistent with the adia- 
batic evolution (spherical or jet) of ly^. 

4.0.6. GRB 971208 

GRB971208 is the longest burst ever detected with 
BATSE. The temporal structure of the burst is a sim- 
ple smooth FRED lasting several thousand seconds. The 
emission is soft, with no emission in channel 4 {E > 300 
keV). The spectral parameters tend to favor fast-cooling, 
but not strongly, as the value of p = 4.06 ± 0.04 is unusu- 
ally high. The value of A = 1.48 ±0.04 is well-determined, 
and very far from the value expected for slow-cooling 
(A = 0.5). Although in apparent contradiction to this. 



the C pD pattern for this event (Figure 7) shows a strong 
resemblance to that of the tail of GRB9890923 in Figure 
5. 

4.0.7. GRB 980301 

GRB980301 shows a low-energy slope consistent with 
fast-cooling but also a value oip — 2.48±0.13, remarkably 
consistent with values of observed afterglows. The change 
in slope across the break energy is slightly higher than that 
expected for slow cooling (but within two-sigma). If the 
spectrum is fast-cooling, then we expect A — 0.74 ± 0.13, 
based on the measured high-energy slope. This value is 
within one-sigma of the value A = 0.69 ± 0.13 that we 
derive from the measured slopes. For radiative evolution 
we expect /? = — 1.84± 0.13, while for adiabatic evolution 
we expect (3 = — 1.36±0.13. However, we measure a much 
steeper value of /? = — 2.50^Q;g5, suggesting an evolution 
inconsistent with the hydrodynamics of a spherical blast 
wave, but consistent with that of a jet. The CCD pattern 
for GRB980301 is shown in Figure 6. Although very simi- 
lar to the model pattern, the observed pattern appears to 
be displaced. 

4.0.8. GRB 981203 

The measured low and high-energy spectral indices for 
this e vent are notably different than those of other bursts 



listed in Table 2. The low-energy spectral index is consis- 



tent with the spectral slope below t'm in the slow-cooling 
mode and below Vc in the fast-cooling regime. Inter- 
estingly, for the fast- cooling regime the high-energy in- 
dex is marginally consistent with the spectral slope for 
Vc < < I'm- The direct implication here is that is 
above the BATSE window and has yet to evolve through. 
Hence the value of p is undetermined. The flux in the 
tail w as too weak to follow the evolution of the spectrum 



our analysis, we have identified a subset of gamma-ray 
bursts that exhibit smooth high-energy (~ 25-300 keV) 
decay emission whose spectral properties are very similar 
to that of fast-cooling synchrotron emission that results 
from a power-law distribution of relativistic electrons ac- 
celerated in a forward external shock. The 25-300 keV time 
histories of the high-energy afterglow candidates arc shown 
in Figure 9. The diversity of the time profiles suggests that 
the GRB time history is not necessarily the distinguishing 
feature of external shock emission in the fireball model. 
The diversity also suggests that the afterglow may be dis- 
connected from the burst emission (e.g., GRB960530), or 
overlap the burst emission (e.g., GRB920622). BeppoSAX 
has demonstrated the existence of both cases: overlap or 
continuation of the afterglow onset with the p rompt burst 
emission [e.g., GRB9 70508 (|Piro ct al 19981) , and more 
recently GRB990510 ( plan et al. 2001] )], and the case in 
which the afterglow begins at a later time, disconnected 
from the prompt GRB [e.g., GRB970228 (|Costa et af 



2000|)]. Our analysis further suggests that in some cases 
(e.g., GRB971208) the early high-energy afterglow may 
actually be the burst emission itself. This situation could 
arise if the energy deposition in the internal shocks is too 
low. 

The spectra of a significant fraction of bursts in our sam- 
ple, however, show inconsistencies with the synchrotron 
model. From a catalog of BATSE GRB spectra, we see 
that the low and high -energy spectral ind ices follow well- 
defined distributions (Preece et al. 2000). For example, 
the distribution of low-energy power-law indices given in 
Figure 7 of Preece et al. (2000) peaks near aiow ^ — 1- 
Roughly 200 of the 5500 spectra in the distribution are 
consistent with the expected value of aiow below t^m, or 
about ^ 4%. If we adopt the hypothesis that GRB spec- 
tra are not synchrotron spectra, then on average ~ 4% of 
the time we expect to measure parameters consistent with 
the synchrotron spectrum purely by chance coincidence. 
This implies that we can expect 0.8 events from Ta- 
ble 2 to have aiow = —1.5 {a = —0.5) purely by chance. 
Clearly, our total of eight candidate events exceeds this 
limit. These events are thus likely sources of synchrotron 
emission. 

Recent studies on electron acceleration models for ultra- 
relativistic shocks predict val ues of the electron index in a 
narrow ran ge 2.0 < p < 2.5 ( Gallant ct al. 1999| , Gallant 
et al. 2000). Our values maintain a significant dispersion 



under the assumption that the observed high-energy after- 
glow is equivalent to the high-energy slope of the fast or 
slow-cooling synchrotron spectrum (i.e., ahigh+1 = — p/2). 
Electron indices in Table 2 tend to be steepr on average 
than p — 2.5. One possible alternative for a high value of 
p (p ~ 3) may be due to a shock generated in a decreas- 
external medium that is the result 



ing density, 

of a massive stellar wind ( phevalier 199g , Chevalier et al 



clearly the temporal decay should be very shallow, unlike 
our measured value of /3 = — l-6i;j]Q oo2- This evolution is 
not consistent with the evolution of a spherical blast wave 
into a constant density medium. 



with my reasonable accuracy. From equation 3 and 4, 1999). On the other hand. Sari (2000) pointed out that 



5. DISCUSSION 



The diverse temporal and spectral properties of GRBs 
leave their origin open to different interpretations. From 



there is no reason why the value of p should be different 
for wind models. Note that several values of p in Table 3 
are below p = 2. Hard electron indices (1 < p < 2) have 
recently been rep orted for the jet model of GRB00301c 
(Panaitescu 2001) assuming a broken power-law electron 
energy distribution. Similarly, a jet interpretation for 
GRB010222 would also require a flatter electron index. 
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Fig. 9. — Time histories (25-300 keV) for the eight high-energy afterglow candidates. 
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Table 4 
Jet Candidates 



GRB 


Trigger 


P 


/3 


910602 
920622 


257 
1663 


2.00 ±0.01(1) 
2.96 ±0.26 


1 74+0.72 
^- '^-0.11 
> 2.98 


920813 
970302 


1807 
6111 


2.02 ±0.02 
1.46 ±0.10 


2 10+''°5 
^•^"-0.04 

1 4Q+0-93 
^•^^-0.11 


970411 
980301 
990316 


6168 
6621 
7475 


2.20 ±0.08 
2.48 ±0.13 
3.04 ±0.09 


2 06+° '^^ 

Z.UD_g 20 

> 2.20 


990322 
990415 


7488 
7520 


0.84 ±0.06 
2.56 ±0.16 


r, 07 + 0.02 

U-«'_o.oi 
> 2.23 


Pre-break Jet Candidates 


910131 


2151 


1.78 ±0.11 


n 71 +0-31 
J^-0.11 


940419b 
960530 


2939 
5478(2) 


3.32 ± 0.55 
2.96 ± 0.22 


> 1.75 
1 40+0-77 


970925 
971208 


6397 
6526 


2.56 ±0.13 
4.06 ±0.04 


^•■^^-0.01 



(1) This value assumes the spectrum is slow-cooling. 

(2) First emisison episode of GRB960530. 



based on analysis of the B eppoS AX daXa, ( In'T Zand et al_ 
20010 However, In'T Zand et al. (2001) show that the 
slowing of the ejecta into the non-rclativistic regime yields 
p — 2.2. From a theory viewpoint, Malkov (1999) has 
shown that it is possible to obtain a hard electron elec- 
tron distribution in Fermi acceleration models. While the 
model predictions for the range of electron indices appear 
somewhat uncertain, the observed dispersion in p values 
may be strongly linked to the accuracy of the fitted value 
of Ohigh- Systematic effects may play a role that intro- 
duces a bias toward steeper values in the estimation of the 



high-c nergy power-law index. 



Our analysis shows that the tail temporal decays are 
well-described by a power-law with a mean index {(3) = 
—2.03 ± 0.51. While the spectral parameters of approx- 
imately 20% of the decays in our sample are in gener- 
ally good agreement with the synchrotron spectrum, the 
temporal evolution, in general, does not agree well with 
the evolution of a spherical blast wave in a homogeneous 
medium. There are alternatives that might explain this 
deviation: first, we only considered fully radiative or fully 
adiabatic evolution. More than likely the fireball is neither 
fully radiative nor fully adiabatic throughout its evolution, 
although at the very early stages the evolution may nearly 
be fully radiative while at latter stages the evolution is 
completely adiabatic. Bottcher and Dermer (2000) have 
considered the early afterglow regime with the intermedi- 
ate cases: partially radiative or partially adiabatic blast 
waves. They find that the temporal decay of the spectral 
flux in the fast-cooling regime is a function of e = eeCrad, 
where trad is the fraction of energy radiated by the accel- 

_1 /, n (! + ') 

erated electrons, or (x v ' t for < v < 

_ ,„ 2(l + ,) + 6(p-l) 

and Fy (X u ^' t for v > i/i^. Higher efficien- 

cies therefore produce steeper temporal slopes. However, 



as can be seen for e — 0.8, we obtain F^ v ^/'^t for 
Vc < V < v^. This is steeper than the expected value of 
— 1/4 from equation 4 but not steep enough to match the 
discrepancies in our observations. 

Another possibility to consider is a jet-like geometry 
or collimated outflow, as opposed to the simple spheri- 
cal blast wave. The break in the light curve may occur 
at early times after the initial shock, as in the case of 
GRB980519, where evidence exists for a break to a steep 
decay that apparently occurred during the few ho urs be- 
tween th e GRB and the flrst afterglow detection ( ^ari et 



al. 199£ ). Rhoads (1999) has shown for adiabatic evolution 



that the time of the break in the observer's frame goes as 
t(, cc 9^. A very early break therefore requires a very small 
Oc- If Oc < Ob initially, then the slope is steep from the 
start. This implies one of two possiblities: (1) very narrow 
emission spots, or "nuggets" , within a narrow coUimation 
angle, or (2) a very small value of F such that F"-'^ > 9c- 
The second option is not likely since the observed emis- 
sion is in the keV to MeV range and oc F^, requiring 
a high Lorentz factor. A list of events from Table 2 and 
3 with comparable electron indices and temporal indices 
(as required for a jet-like blast-wave) are presented in Ta- 
ble 4. Note that although the value of p for GRB990322 
is low, the temporal decay does follow the t^'P relation. 
Events with values of (3 shallower than p may be events 
in which the break occurs at some later time after the 
GRB. We categorize these events as pre-break jet candi- 
dates. Three of these events (GRB940419b, GRB960530, 
and GRB970925) have some spectral properties charac- 
teristic of the synchrotron spectrum (see last column of 
Table 2). The flve bursts labeled in Figure 3 are also candi- 
dates for jet outflows. More importantly, note that three of 
these events (GRB920622, GRB980301, and GRB970411) 
are strong candidates because they belong to the group of 
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high-energy afterglow candidates presented in §4 that were 
selected based on their spectral properties. The spectra of 
the r emaining two events (GRB920801 and GRB931223) 



cay emission that softens relative to the initial burst emis- 
sion, suggesting the existence of early high-energy after- 
glow. Other studi es using PHEBUS (Tkachcnko et al 



are oi ly marginally consistent with synchrotron emission 
from an external shock. 

6. SUMMARY AND CONCLUSION 

Our temporal and spectral analysis of the smooth ex- 
tended gamma-ray decay emission in GRBs has shown ev- 
idence of signatures for early high-energy afterglow emis- 
sion in gamma-ray bursts. The extended decay emission is 
best described with a power-law function F^, oc rather 
than an exponential, similar to the results of Ryde and 
Svensson (2001) who studied the decay phase of a sample 
of GRB pulses with a broad range of durations. From our 
sample of 40 events, we find {(3) « —2 for long, smooth 
decays. Color-color diagrams have provided a qualitative 
interpretation of the burst spectral evolution and allow a 
simple comparison with the evolution expected from the 
synchrotron model as well as comparison of spectral evo- 
lution among GRBs. The CCD patterns and the spectral 
analysis indicate that ~ 20% of the events in our sample 
are consistent with synchrotron emission expected from an 
external shock. Interestingly, three of these events have de- 
cay rates consistent with that expected from the evolution 
of a jet, Fi, ~ t~P. Because the break is essentially at the 
onset of deceleration, the jet must, at least, be very nar- 
row, since 6c < l/F. Table 4 suggests that in some cases 
the break occurs at a later time, so that the prompt emis- 
sion we observe is pre- break, 9c > l/T and consistent with 
spherical geometry. A possible scenario is one in which the 
ejecta is very grainy, where the nuggets in the ejecta are 
smaller than l/T, similar to the model discussed by Heinz 
and Begelman (1999). Huang et al. (1999) (see also Wei 
and Lu 2000) have shown that the break in the light curve 
is more of a smooth transition due to the off-axis emission 
of a jet with no angular dependence. The steep light curve 
can only occur if the angular size of the nugget is less than 

1/r. 

Connaughton (2001) has investigated the average late- 
time temporal properties of GRBs observed with BATSE 
and found statistically significant late time power-law de- 



2000D and APEX ( [Litvinc et al. 200C| ) bursts show sim- 



ilar behaviors in late-time GRB light curves. Collectively, 
these studies strongly suggest that the afterglow emission 
may overlap or be connected to the prompt, variable burst 
emission. On the other hand, it is clear that not all GRBs 
exhibit such behavior. In some cases, the initial gamma- 
ray flux from the external shock may simply be too low to 
detect (e.g., see Figure 4 in Gibhn et al. 2000). In other 
cases, the bulk Lorentz factor may be too low to generate 
the gamma-ray photons upon impact with the surrounding 
medium. 

As the number of afterglow/counterpart detections in- 
creases, the relationship of the afterglow emission to the 
gamma-rays released in the initial phase of the burst 
can be studied systematically. The capabilities of Swift 
(Gehrels 2002) will allow broad spectral coverage using 
three co-aligned instruments (BAT, XRT, and UVOT) 
during the gamma-ray phase and early afterglow phase of 
the burst and facilitate the distinction between the GRB 
and the onset of the afterglow based on temporal and 
spectral information. With well-constrained spectral and 
temporal parameters in hand, plots of temporal index vs. 
spectral index can be readily constructed and thus provide 
information on the geometry of the fireball and definitively 
test the internal/external shock model for GRBs. 
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